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We study the phase structure of the Nf = 2 flavors QCD with imaginary chemical po- 
tential with the use of the clover-improved Wilson quark action and renormalization-group 
improved gauge action. We calculate the Polyakov loop on 8 3 x 4 lattice for (jS , ft") = 
(1.8,0.1411), (1.9,0.1388), (1.95,0.1377) with ju 7 = 0.2618. We find that the phase of the 
Polyakov loop shows a two-state signal indicating the first order phase transition. This transi- 
tion occurs in the vicinity of j3 = 1.9, which corresponds to the temperature T/T pc = 1.08. 
We also present a reduction formula for the quark determinant of the Wilson fermion. We discuss 
the feature of the matrix reduction formula. 
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1. Introduction 

It is of prime interest to understand QCD at finite temperature and density. The lattice QCD 
with a non-zero baryon density suffers from a sign problem, where the determinant of the quark 
matrix becomes complex and Monte Carlo simulations become unavailable. Because of this, it is 
still challenge to understand the QCD at finite density. Several approaches have been investigated 
to overcome or circumvent the sign problem: multi-parameter reweighting, canonical approaches, 
imaginary chemical potentials, and so on. 

If the chemical potential is pure imaginary }i = pL R + iHi(HR = 0), the quark determinant be- 
comes real due to the identity A(/i/) t = Y5^(^i)Y5- Thus, there is no sign problem in this case, and 
the standard Monte Carlo algorithm can be applied to. The imaginary chemical potential was in- 
vestigated with the use of the staggered (KS) fermion by De Forcrand and Philipsen [|l|] and D'Elia 
and Lombardo [Q, 0]. Wu, Luo and Chen considered the imaginary chemical potential using the 
standard Wilson quark action 

The study of the phase structure in (j8,/^/)-plane provides us with an understanding of the 
phase structure in (j3,£i#) -plane through analytic continuation. In the study of the imaginary chem- 
ical potential, it is important to consider the Roberge-Weiss periodicity In the presence of the 
quark, the Z(3) symmetry is explicitly broken. However, if the chemical potential is pure imagi- 
nary, the Z(3) symmetry is maintained via a translation of pLj as 

where k is an integer. Roberge and Weiss showed the existence of a phase transition at high tem- 
perature, where the phase of the Polyakov loop is an order parameter, while this phase transition 
does not occur at low temperature. A critical line corresponding to this phase transition is ^ = 
with a critical endpoint Tg, see Fig |l[ 




7T/3 2ti/3 \ij/T 

Figure 1: A schematic figure of the phase structure in (T,ji]/T) plane. 
Vertical solid lines show the critical lines of the Roberge-Weiss phase 
transition. Vertical dashed lines shows the Roberge-Weiss periodicity. 
The crossed symbols corresponds to points where the present simula- 
tions are performed. 

In this study, we present two results. First, we report a preliminary results on the study of 
the phase structure of the Nf = 2 flavors QCD with an imaginary chemical potential by using the 
clover-improved Wilson quark action and renormalization-group improved gauge action on 8 3 x 4 
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lattice for (j3,K) = (1.8,0.1411), (1.9,0.1388), and (1.95,0.1377) with the imaginary chemical 
potential = 0.2618. Second, we present a reduction formula for the quark determinant of the 
Wilson action, which is an extension of the reduction formula given in Ref ^, §]. 

2. Phase Structure with Imaginary Chemical Potential 

2.1 Formulation 

We employ the renormalization group improved gauge action []|] 

S s = § x 1 loop)+ci£(l x 2 loop)] , (2.1) 

with c\ = —0.331 and cq = 1 — 8c i and clover-improved Wilson action, where the quark matrix is 
written as 



— fC 



— S.xjCswK E a iivF^ v - (2.2) 

ft<V 

For the coefficient of the clover term C$w> we us e a result obtained in the one-loop perturbation 
theory [fill], C sw = (1 — 0.8412j3 J ) 3 / 4 . We calculate the Polyakov loop (L) = |L|exp(/0), where 



the phase is an order parameter of the Roberge-Weiss phase transition [|p. 
2.2 Results 

We perform simulations on 8 3 x 4 lattice by using the standard hybrid Monte Carlo algorithm. 
In the present study, we carry out the simulations for three set of parameters (jS , fc) = (1.8,0.141 1), 
(1.9,0.1388), and (1.95,0.1377). All for these three cases, we use a fixed value of the imaginary 
chemical potential jij = 0.2618, which is on the line of the Roberge-Weiss phase transition Hi/T = 
7i/N c . These parameters correspond to lines of the constant physics with mps/my = 0.8 in the 
absence of the chemical potential [p||. We set a step size of the molecular dynamics to be 8z = 0.02 
and the number of the molecular dynamics to be 50, which gives the length of a molecular dynamics 
trajectory to be one. We generate 10, 000 trajectories and measure the Polyakov loop for each 
trajectory. 

Figures || show the history of the Monte Carlo simulations for the absolute value and the phase 
of the Polyakov loop for /3 = 1.9. The phase of the Polyakov loop <p shows the two-state structure 
with 0=0 and —n/3. This indicates the first order Roberge-Weiss phase transition. 

Figure || shows the distribution of the phase of the Polyakov loop 0, where the first 5,000 
steps are removed as thermalization. Following Ref. [fn]], three sets of the parameters correspond 
to the temperature T/T pc = 0.93 for (]3,k) = (1.8,0.1411), T/T pc = 1.08 for (1.9,0.1388) and 
T /T pc = 1.20 for (1.95,0.1377), where T pc is a pseudocritical temperature. It follows from the 
figures that at low temperature (/3 = 1.8), the system has one vacuum with the value of the order 
parameter (j) ~ 0.3. As temperature becomes higher (j8 = 1.9), two-state behavior appears with the 
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Figure 2: Monte Carlo History of the Polyakov loop for (P,K,jj,j) = 
(1.9,0.1388,2.618). The left panel shows the absolute value \P\ and the right 
panel shows the phase 0. 

Distribution of ()) 




Figure 3: A histogram of the distribution of the phase of the Polyakov loop (j). 
The dashed (black), solid (red) and dotted (blue) lines correspond to (j3 , fc) = 
(1.9,0.1411), (1.9,0.1388) and (1.95,0.1377), respectively. We use the same 
value of the imaginary chemical potential jj.] = 0.2618 for these three sets of 
the parameters. 



two vacua 0=0 and <p = —7i/3, although the peak structure is not sharp due to the small lattice 
size. As temperature becomes much higher (j8 = 1.95), the vacuum shift to <p = 0. We find the 
onset of the Roberge-Weiss phase transition is in the vicinity of j8 = 1.9 for which T /T pc = 1.08. 
This is almost consistent with the previous result obtained by the standard Wilson quark action m\. 



14, 15] reported larger values Tg 



On the other hand, the phenomenological approaches [ ]12| , [H 
200 ~ 210 MeV. This dicrepancy is probably caused by the fact that we employed large quark mass 
nip S /my ~ 0.8. 



3. Gibbs formula for Wilson fermions 



In this section, we derive a formula which is indispensable if we use Wilson fermions for 



4 



Wilson fermions with imaginary chemical potential 



Keitaro Nagata 



the multiparameter reweighting by Fodor and Katz (i.e., using no Taylor expansion) , and also 
for the canonical expression. Danzer, Gattringer and Liptak studied a decomposition of the quark 



determinant for Wilson action[[l^, |17p. Here we derive an alternative expression. 
We can write Wilson fermion matrix as 



A = B — z~ l K(r — y A )V - zx(r + y A )V\ 



(3.1) 



where 



B(x,x') = 8 x , x ,-K^Ur-Ym{x)^, x+ i + {r + Yi)Ul(x>)8^_~X 
+ {Clover), 



and 



V(x,x>)^U 4 (x)8 x , x+i . 
We introduce z as z = . Now we rewrite A: 

detA = Aet{B-zr l K(r-y A )V - z*{r + y A )V r ) , 
-EV-z{-K(r+y A )) I 
K(r-y 4 )V 2 -z 



(3.2) 
(3.3) 



/detV 



-N 



-BV I 

K(r-y 4 )V 2 



■K{r + y A ) 
/ 



(3.4) 



Here Af is a rank of the block matrices, such as B and V, N = N c x 4 x N x N y N z N t . By exchanging 
the columns and raws, this matrix now reads 



detA = z~ /v det(r-zS). 
Here we describe the matrices T and S as block matrices in time-plane, 



(3.5) 



(o 


h 


•■• 











h ■■■ 











••• 








• • • tN, -2 











... o 




Vn, 





... o 






, s 











where 



-BiVu+i 1 



Kir-y^Vt-uVu+i J 
Each ti and s is (4N c N x N y N z ) x (4N c N x N y N z ) matrix. 



-K(r + y A ) 
/ 



(3.6) 
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Even though the matrix S does not have an inverse, the formula (3.5) can be evaluated by the 
generalized eigen value problem, 

TX = zSX. (3.7) 

Namely the generalized Schur decomposition^] tells us that there exist unitary Q and Z such that 
Q^SZ and Q^TZ are upper triangular. We write their diagonal elements, a, and j8,, respectively. 
Then 

det(r-zS) =det(<2Z t )[J(a,--zA-)- (3.8) 

This is a Gibbs formula for Wilson fermions. However, we evaluated a, and j8, using LAPACK 
routine, and there is an accuracy problem. Then we go further to obtain a satisfactory formula. 
Let us rewrite the determinant of T — zS as 



det(r-zS) 








... o 


-zs 


—zs 


h 


■•• 








-zs 


h ■•• 











—zs ■■■ 








■ ■ ■ 











■■■ -ZS 





We multiply matrices Q\ and Q2 from the right 



Qi 
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/ 
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••• 
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••• 











/ ••• 















/ ••• 



































1 ) 




V 










1 J 



Here detgi = det<22 = 1- Then, we obtain 

det(r-aS) = det(T-z.S)Q i Q 2 , 

= \tN,\ 



X 



h 










-zs 


h 




-zst^ l zst N ^ zs 





-zs 


■ ■ ■ fy,-2 










••• — zs 


fiVf — 1 



\t Nt \ X |^| X 



h 




-zst 1 zst Nt zs 


-zs 
















■■■ -ZS 


tN,-l 



(3.9) 



(3.10) 



(3.11) 
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We go further recursively. 

det(r-zS) = det(r-zS)gig 2 ---ftv,-2, 



\t N , \ X I?!] X •■• |^_ 3 1 
tN,-2 -ZSt^ZS-'-t^ 1 ZSt Ni l ZS 

-zs t Nt -\ 



(3.12) 



Since 



tN,-2 ~ZSt Nt l _ 3 ZS---t x l ZSt Ni l zs\ _ 
-ZS t Nf -\ J 

t N ,- 2 o\f I -t^ t 1 _ 2 zst^_ 3 zs---t{ 1 zst^ 1 zs 
Oil \-zs t Nt -i 



(3.13) 



we obtain 



det(T-aS) = |^| x \h\ x ■•■|%_ 2 | 

X |?iV ( -l -Z5'? iVf 1 _ 2 Z5^ 1 _ 3 Z5'---?lZ5'^ 1 Z5'|, 
= X |*i | X ■•■|*JV,_ 2 | X l%-l| 

x \I - t N ^_ { zst N ^_ 2 zst Nt l _ 3 zs---h zst Nt l zs\, 

= \P\ x \1 - z N 't N ^ l st N ^ 2 st N ^s- ■ -hst^sl (3.14) 

where P = t\t2 ■ 

Therefore, the ratio of the fermion determinant with z = to that with z = 1, i.e., \l = is 

detA( M ) __ N det (I -z N 'Q) _ ^ 



detA(/i=0) det(/-<2) 
Here Q is a matrix of L x L with L = 2 x 4 x N c N x N y N z and is given as 



2 = fy-l st N,-2 St N,-3 s " • " Z l sZ A', 1 S " ( 3 - 16 ) 



If we can diagonalize 2 — > diagjgi , q 2 , ■ ■ ■ , then 



detA( M ) .-^ niiCl-A?/) 

~ z ttL n — ~T' (3 - 17) 



detA(jU=0) ltf =1 (l-</z) 
Although it contains inverse matrix calculations, t~ l , the matrix 2 does not depend on 7V f . 

4. Summary 

We have investigated the phase structure of the QCD with the imaginary chemical poten- 
tials. We employed Nf = 2 flavors QCD with the renormalization-group improved gauge action 
and the clover-improved Wilson quark action. We performed simulations on 8 3 x 4 lattice for 
(P,k) = (1.8,0.1411), (1.9,0.1388), (1.95,0.1377) with = 0.2618. We found that the phase of 
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the Polyakov loop showed the two-state signals at (j8, k) = (1.9,0.1388), which indicates the first 
order phase transition. This phase transition occurs in the vicinity of /3 = 1.9, which corresponds 
to T /T pc = 1.08. In order to understand the phase structure of the QCD, such as the pseudocritical 
line J3 C (;U) and the endpoint of the Roberge-Weiss phase transition, the simulation for other values 
of (j3, K, jif) is in progress. 

AN is supported by the Grant for Scientific Research [(C) No.20340055] from the Ministry of 
Education, Culture, Science and Technology, Japan. 
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